Periodic vacuum and particles in two dimensions 



o 

Q 

in 



X 



Marianne Dufour Fournier a *, Janos Polonyi 06 ''' 
""Laboratory of Theoretical Physics, Louis Pasteur University 
3 rue de I'Universite 67087 Strasbourg, Cedex, France 
b Department of Atomic Physics, L. Eotvos University 
Puskin u. 5-7 1088 Budapest, Hungary 
(February 1, 2008) 



ON 
ON 

ON ■ Abstract 



Different dynamical symmetry breaking patterns are explored for the two 
dimensional cf) 4 model with higher order derivative terms. The one-loop sad- 
dle point expansion predicts a rather involved phase structure and a new 
Gaussian critical line. This vacuum structure is corroborated by the Monte 
■ Carlo method, as well. Analogies with the structure of solids, the density 

CN ■ wave phases and the effects of the quenched impurities are mentioned. The 

, unitarity of the time evolution operator in real time is established by means 

<^> of the reflection positivity. 
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I. INTRODUCTION 



Ph. The condensates occuring in Quantum Field Theory are usually homogeneous and are 

composed of particles with vanishing momentum. In this manner the momentum of the 
excitations is conserved even when a particle is borrowed from or lended to the vacuum. 
But the momentum conservation observed experimentally at finite energies is actually not 
incompatible with certain inhomogeneous vacuua so long as the momentum of the con- 
densed particles is beyond the observational range. The elementary excitations in solids are 
described by the Bloch waves which can be rearranged into different sub-Brillouin zones in 
such a manner that the Bloch momentum, the momentum counted from the center of the 
sub-Brillouin zone, is conserved. More formally, the presence of a crystalline ground state 
restricts the translations as symmetries such that the primitive unit cells are mapped into 
each other. The conserved quantum number due to such a restricted symmetry group is 
the Bloch momentum. The umklapp processes which take place at the length scale of the 
primitive unit cell change the sub-Brilllouin zone and can be interpreted as a change of 
the type of the excitations. Returning now to Quantum Field Theory, one might send the 
size of the primitive unit cell to zero. If this is possible then the space-time structure of 
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the momentum non-conserving umklapp process is not resolved by finite measurements and 
their interpretation as a flavor changing process becomes compatible with the experiments. 

In order to gain more insight into the role the inhomogeneity of the vacuum plays in 
forming the dynamics of the excitations we consider a generic model with higher order 
derivative terms in the action for a scalar field in two dimensions and present its phase 
structure, the excitation spectrum and the particle content when the vacuum possesses a 
modulation and becomes inhomogeneous. We use the saddle point approximation in the one- 
loop order in the analytical computation and the Monte Carlo for the numerical simulation to 
have a more complete picture. The higher order derivatives lead in general to the appearance 
of additional particles with negative norm and complex energy. It is pointed out that our 
model possesses the reflection positivity which in turn assures the existence of the positive 
norm Hilbert space and the unitarity of the time evolution operator in Minkowski space-time. 

The organization of the paper is the following. Section II contains our motivation in 
choosing the model investigated. The tree level vacuum is identified in Section III. The 
action is rewritten in terms of the Bloch waves corresponding to the different periodic vacuua 
in Section IV. We diagonalize the quadratic part of the action and determine the elementary 
excitations for the simplest inhomogeneous vacuum in Section V. Section VI contains the 
demonstration of the one-loop renormalizability of our model. The analytical results are 
compared with a Monte Carlo computation in Section VII. The issue of the unitarity is 
discussed in Section VIII. Finally Section IX is for the conclusions. 

II. THE MODEL 

Our model is an extension of the Landau-Ginzburg model for a scalar order parameter 
by adding higher order derivatives to the action, 

S[<f>(x)} = J d- x{^(x)/c(^n)^(x) + V((f>(x))}, (1) 
where the kinetic energy contains the functions 

K{z) = 1 + c 2 z + c 4 z 2 , y(0) = ^V + ^> 4 , (2) 

and A is the ultraviolet cutoff. 

As far as the dimensionless parameters c 2 and C4 are concerned, we have two different 
motivations for their introduction. One is based on the fact that we are always confronted in 
Nature with effective theories where the high energy particle exchanges generate a number 
of operators in the action which are perturbatively non-renormalizable. As an example, 
consider a renormalizable model for a heavy and light particle, described by the fields $ and 
0, respectively and the bare action S [(f), <&]. The effective action for the asymptotic states 
below the threshold of the heavy particle is given by S e ff[<f>] = S [(j), $ = 0] + AS[<fi], where 
AS[<f>] contains all effective vertices generated by the heavy particle exchange processes. In 
this manner we can never be sure that the action corresponding to the interactions in a 
given energy range is actually limited by the renormalizability even thought the "Theory 
of Everything" is supposed to be finite or renormalizable. The higher order derivative 
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terms of our model may arise from AS*[0]. The decoupling theorem, helps us out from 
the problem of a too general action with nonrenormalizable terms by asserting that the 
nonrenormalizable coupling constants are small, being suppressed by the power of the light 
and the heavy particle mass ratio. The model ([!]), a slight extension of an effective theory 
for the Higgs boson || , retains some of the suppressed nonrenormalizable terms arising from 
a hitherto unknown super heavy particle exchange. 

The question, left open by the conclusion of the decoupling theorem, and which motivates 
the present work is whether the smallness of the non-renormalizable coupling constants is 
really sufficient to render them unimportant in the effective theory. We shall find that certain 
higher order derivative terms may become relevant when their coupling constants exceed a 
threshold value. In other words, we might be forced to consider non-renormalizable terms 
in our effective theories if the heavy particle is not exceedingly far from the observational 
energy. The other motivation to study the model (jl]) with non-renormalizable terms is the 
suspicion that the perturbatively non-renormalizable terms might turn out to be relevant 
by a non-perturbative mechanism and actually allow the removal of the cutoff. If this is 
happened to be the case then our model with non-vanishing C2 and C4 is as justified as the 
usual one whose action is quadratic in the gradient, a possibility motivating the authors of 
ref. §. 

What kind of heavy particle exchange is behind the higher order derivative terms of our 
action? As mentioned above, our interest is in theories with inhomogeneous vacuum. Such 
a ground state which is modulated with a period length X vac is the result of a force which 
is attractive for x > X vac and repulsive when x < X vac . When the particle whose exchange 
generates this force is eliminated, its effects are kept in the choice of the vertices in AS[<f>]. 
The momentum independent ultra-local (i.e. non-derivative) terms which contribute to the 
local potential in the action can not generate such a strong distance dependence in the 
interaction. But it is easy to see that the higher order terms in the derivative are just for 
this role, to lower the action for modes whose characteristic momentum scale is p ~ 1/X vac . 
In fact, consider the eigenvalues of the second functional derivative of the action evaluated 
in the trivial vacuum < <p{x) >= 0, 

e (p) = ^(0)+^(-^ 2 ). (3) 

The 0{p A ) term produces a non-trivial local minimum at p — p m in ~ ^vac f° r c 2; c 4 > 0. 
Thus C2 and C4 correspond to a van der Waals force. We shall go further in this work and ask 
what happens when c<i reaches so large values that the minimum of the dispersion relation 
e(p) turns out to be negative. The corresponding vacuum will be the condensate of particles 
with momentum around p min and will display the period length X vac ~ p^L- 

The hand-waving argument to retain the higher order derivative terms of ([!]) from the 
multitude of other contributions in A5[0] is the following. Let us start with C4 = 0, when 
the vacuum is homogeneous for C2 < 0. On the contrary, for C2 > an instability opens 
by increasing the momentum of the condensed particles and the ultraviolet cutoff stabilizes 
the vacuum || where particles with momentum at the cutoff are found. Models with such 
an instability were studied in refs. ||, [BJ and || in three and four dimensions. If the 
inhomogeneous vacuum is supposed to be formed at momentum scales below the cutoff 
then we need another stabilization mechanism. For this end we retain the 0(<9 6 ) term with 
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C4 > 0. We believe that the higher order terms in the gradient will modify the shape of the 
saddle point only leaving the qualitative features of the inhomogeneous vacuum unchanged. 
In other words, the kinetic energy is generic. The present work can be considered as the 
continuation of Ref. |J where a new gaussian ultraviolet fixed point was found in the one- 
loop approximation for C4 = 0. We simplify in this work the issue of the renormalizability 
by choosing lower dimension, d = 2, but the multitude of different phases is explored by 
allowing C4 7^ 0. 

The allusion made above at the Landau- Ginzburg model is based on the similarity of 
the functions V(u) and K(u). The nontrivial absolute minima provide the mechanism of the 
spontaneous (p = 0) or dynamical (p 7^ 0) symmetry breaking generated by the potential 
in the internal space or the kinetic energy in the external (and internal) space, respectively. 
Our interest in this work is to explore the different dynamical symmetry breaking patterns 
provided by the generic action (D) and to suggest a mean-field treatment for the phase 
transitions with modulated ground state JQ. 

The questions addressed here and their tentative answers have certain relevance both in 
Solid State and High Energy Physics. The nontrivial, periodic vacuum generated by the 
higher order derivative terms may offer a new point of view in understanding the origin of 
the crystalline structure in solids. In fact, consider the coupled system of electrons, ions and 
photons. It is already an effective theory because the lower lying electrons of the ions are 
represented by the insertion of different charge distributions and form factors for the ions. 
We introduce chemical potentials for the electrons and the ions in order to realize electrically 
neutral matter with finite density. Finally we eliminate the heaviest degrees of freedom, the 
ions. This is opposite to the usual Born-Oppenheimer approximation but it generates a 
local effective interaction. This effective theory for electrons and photons contains higher 
order derivative terms for the photon field which yields the periodic, crystalline ground state. 
Another appearances of this mechanism where strong van der Waals forces are acting are the 
antiferromagnets and the charge or density wave phases. In the latter the effective theory 
is obtained by eliminating the valence electrons and has periodic tree level ground state || . 
The massless case with e{p min ) > may as well represent a superfluid system with p min as 
the roton momentum ||. In this context our quantitative argument showing the relation 
between the van der Waals forces and the higher order derivatives in the action represents a 



simple, effective theory motivated alternative of Feynman's argument about the rotons |L0 



The instability leading to the formation of an inhomogeneous ground state when e(p m i n ) < 
is a quantum phase transition where a non-classical soft mode [[□]] shows up at p ~ p m in- 

The periodic vacuum of our model supports frustrations for certain choice of the coupling 
constants. To see this we write the lattice regularized version of ([TJ) in the form 

s[<K*)] = £{-^W D )<K*) + vMx))} 

= El^O) \a<k*) + E ( B ^ x + e ») + <f>ip - ^)) 

+ C(<P(x + 2e M ) + <f>(x - 2e M )) + D{<f){x + 3e M ) + <f>(x - 3e„))) 

+ E ( E (4>(x + e M + e » ) + <t>{ x + e M - e„) + (j)(x + e v - e M ) + <j>{x - e M - e„)) 
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+ F(<f>(x + 2e M + e„) + 4>(x + 2e^ - e u ) + 4>(x - 2e M + e u ) 

+ 4>{x - 2e M - e„))) (4) 
+ G ^ (c/>(x + e M + e v + e p ) + 3c/>(x + e M + e„ - e p ) 

+ 3c/>(a; + e^-e v - e p ) + 0(x - e M - e y - e p )^ + ^0 4 (x) J, 

where A = m 2 /2 + d - (2d 2 + d)c 2 + (4c/ 3 + 6c/ 2 )c 4 , 5 = -1/2 + 2c/c 2 - (6d 2 + 3d - 
3/2)c 4 , C = £/2 = -C2/2 + 3cfc 4 and D = F/3 = G = -c 4 /2. The possibility of having 
either sign for these coefficients indicates the competition between the nearest- and beyond 
nearest-neighbor interactions and the possible presence of frustrations as lattice defects of 
the periodic vacuum. When the amplitude of the periodic vacuum is large the motion of the 
frustration is rather slow and the model offers a semiclassical description of the quenched 
disorder. 

The usual strategy in High Energy Physics is to introduce a field variable for each 
particles. But models with inhomogeneous vacuua may display more involved particle- 
field assignments by exploiting the non-trivial dynamics at the ultraviolet cutoff scale. In 
fact, there are several dispersion relations and particle like-elementary excitations in Solid 
State Physics, such as the acoustical and the optical phonons and the massless or massive 
excitations in antiferromagnets. In both examples the dynamics is rather non-trivial at the 
ultraviolet cutoff. Can we keep such a more unified description of several particles by means 
of a single quantum field in a renormalizable model? In this case the dynamics at the cutoff 
can be pushed at infinitely high energies and it is not obvious that the construction converges. 
Our answer to this question is affirmative up to the one-loop order of the perturbation 
expansion. Another interesting aspect of the model considered is the possibility of breaking 
continuous external symmetries without generating massless Goldstone modes. In fact, 
suppose that we use say lattice regulator which breaks the continuous translation invariance 
and the momentum of the condensed particles is close to the cutoff, in which case the 
continuum description is not applicable at the length scale of the vacuum and there is 
no reason to expect gapless excitations. This is another manifestation of the apparent 
homogeneity of the modulated vacuum with shrinking period length. 

There are several works devoted to the applications of models with higher order deriva- 
tives. The canonical formalism for Lagrangians with higher order derivatives has been 
worked out in ref. 0] and been applied to QED fl3f , [fL4}| , gravity |16| , the Higgs sector of 
the Standard Model and Chern-Simons theories [17]. The present paper is close in spirit 
and the method followed to the nonperturbative lattice study of the effects of the higher 
derivative terms for a scalar particle ||. But contrary to that paper where the continuum 
limit is sought in the symmetrical and the ferromagnetic phase by keeping both the particle 
mass and the scale induced by the higher order derivative terms finite we shall consider the 
phase with inhomogeneous vacuum whose length scale does shrink to zero with the cutoff. 

The higher order derivative terms may generate new tricritical point, called the Lifshitz 
point. One can see this in the Euclidean field theory given by the lagrangian 



!(^) 2 +fc< 



) 4 + f <p 2 + 10 4 + |0 e . 



(5) 
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It is well known that the model C2 = 0, ge > possesses a tricritical point which is at 
92 — 94 — in the mean-field approximation. For g A > the phase transition following 
the sign of gi is changed is of second order. When g A < then the crossing of the line 
92 = 5gl/8g e corresponds to a first order phase transition whose latent heat vanishes as 
(74 — ^ _ . A similar tricritical point is generated by the competing derivative terms for 



ge = at gi = C\ = 0, according to an e-expansion like renormalization group method [18 



and the solution of the spherical model [19]. Phase transitions of different type coincide 
here. The change of the sign of g<i induces a second order phase transition when c\ > 0. 
(4>(x)) is inhomogeneous for c\ < 0, and its wave vector approaches zero as C\ — > 0~. All 
of these studies deal with the quantum or statistical physics of the fluctuations around the 
trivial vacuum, (<p(x)) = 0. We go beyond this level by allowing the condensation in the 
ground state and considering the fluctuations around a periodically modulated mean field 
(0(x)>. 

The interference of the cutoff and the derivative terms leads to commensurate- 
incommensurate transitions | 20|j , as well. The simplest 1+1 dimensional model of these 



transitions is based on the lattice potential energy 



+ 1) _ + 5 f - 9 - cos 27r0(a;) 



2 V " V ' " v ' ' 2 b 



(6) 



For large g the value of the field tends to take integer multiples of b, the period length of the 
second term in the potential. When z is chosen to be large then the average increment of the 
field at neighboring sites, 5, is close to 5. The configurations which minimize the energy are 
governed by the competition between the conflicting requirements of the two contributions 
in (]|). The ratio S/b plotted against 5/b exhibits a characteristic, discontinuous structure, 
called the devil's staircase. Similar competition is expected between the period length of 
the vacum X vac introduced at eq. (§) and the lattice spacing in the solution of the equation 
of motion of (f|). Our model offers a translation invariant realization of the commensurate- 
incommensurate transitions. Furthermore the sign of d/j,{4>(x)) is given by the dynamical 
breakdown of the space inversion symmetry instead by the explicit symmetry breaking term, 
6, in®. 

We close this section by mentioning a formal problem which poses a serious threat in using 
actions with higher derivative terms ||13|| . The inverse propagator (|3]) has several roots as the 
function of p 2 . Each of them yields a (not necessarily) simple pole for the energy integrals 
which is passed by the integration contour during the Wick rotation of the perturbation 
series into Minkowski space-time. When the real part of the Euclidean energy at a simple 
pole is non-vanishing then the Wick rotation produces a particle mode with complex energy 
and with amplitude which growths or decreases exponentially in time. Even if the Euclidean 
theory is stable we have to make sure that such runaway amplitudes are canceled after the 



Wick rotation [14]], |]I5 |. Since the original theory, the source of our effective model was 
supposed to have stable vacuum in the Minkowski space-time this problem may only arise 
from the inconsistent truncation of the effective lagrangian ([[]). 

One encounters a similar problem in the Schwinger- Dyson resummation procedure, as 
well. The effective action which approximates the Wilsonian renormalized action in the 
infrared limit, as the cutoff tends to zero, has a non-trivial momentum dependent piece 
which is quadratic in the field, the self energy. When the self energy is considered in a 
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non-perturbative manner then its momentum dependence may create new poles. These are 
usually ignored since they imply that the self energy, coming from the radiative corrections 
is as important as the original p 2 kinetic energy term. But their proper treatment leads to 
the question of the consistence of the truncation of the effective action. 

This point of view opens a general question. All theory we know in Physics is effective 
only and contains higher order derivatives, whatever small coefficients they may have. How 
can we make sure that the studies restricted to the models with the usual kinetic energy 
reflect what happens in the effective theories? The answer to this question involves the 
verification of the usual universality scenario in the mixed framework of the static and 
the dynamical renormalization group, where the fixed point is in the ultraviolet for the 
spatial momenta and in the infrared for the energy. Instead of embarking this important 
but extremely involved issue we restrict ourselves in this work to show that our model is 
consistent in real time and the runaway modes can be excluded in a nonperturbative manner. 



III. TREE LEVEL VACUUM 



We start the semiclassical solution of our model by determining the minimum of the 
lattice action. This is in principle a rather involved numerical problem, the minimization of 
(B). To circumvent this complication we shall seek the tree level, mean-field vacuum in the 
form 



WW = 4>H + 0/i/COS 




M 

K = 2vr— 
N 



(7) 



where the amplitudes 4>h, 4>ih, the relative primes N, 1 < M < N/2, the phase 9 and the 
number of the antiferro magnetic directions, cLaf = 1,2 serve as the variational parameters 
to minimize the action. We have naturally to confine our study into regions far from the 
critical points, i.e. above the Ginzburg temperature in order to apply this method. The 
phase is called para-, ferro-, antiferro- and ferri-magnetic for (p H = = (()jh, 4>h 7^ = 4>ih, 
4>h = 7^ (piH and 4>h ^ ^ 4>jh- The action density, s(cIaf, M, N) = S/L 2 , on a lattice 
L x L with m 2 < is Sq = — m 4 /4A, for d^F — 0. 

For the computation of the action for dAF > we restrict ourselves to the case m 2 < 0. 
The mean-field vacuum configuration is an eigenvector of the lattice box operator, 



a 4>MF(x) = — AdAF sin 2 



K 



)MF\X) - <pH) 



what gives 



UK{U) ( j )MF (x) = M 2 (M, N, d AF , c 2 , c A ){(t)MF{x) - <p H ) 



with 



M 2 = Ad A Fsm 2 



K 



1 - id K.-K111- ( c 2 + 16d^ F sin 4 ) c j 



(8) 



(9) 



(10) 



The dependence of the eigenvalue M. in the parameters M, N, dAF, c 2 , c 4 will be suppressed 
in the expressions below. The action density to minimize is 
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m 
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7 



+ 



C 2 (N) 



- — \m 2 + M 2 

3C 3 (N) 



, „ 3C 2 (N) j2 j2 CAN) . 



(11) 



where we introduced the notation 

C n (N) = 



(12) 



Notice the M-independence of the sum for the relative primes M, N. Direct computation 
gives 



C 2 (N) 
C 3 (N) 

C A (N) 



0= 



cos 20 \ 
2 J 



iV = 2, 
iV > 2, 



2~cUf 

^" AF iV = 3, 
N^3, 







4 / 



(13) 



' ^ 3+4 cos 20+cos 46 ^ dAF N = 2 

N = 4, 

iV = 3 or N > 4. 



3+cos 40 \ 
8 >/ 

I (ff 



<i>lF 



The dependence of C n in A" will not be shown explicitly below. One can see from the above 
construction that the mean-field solution of a model <p e requires the coefficients C n with 
n — 1, • • • ,£. So the limit A" — > oo, the regular dependence on A" sets on for N > £. 
The action density corresponding to different choices of A" is obtained as follows: 
N = 2 : The phase parameter 9 is redundant in this case. We choose 9 = 0, what sets 
C 2 — C4 = 1, and write 



A ■ X + -X ■ B ■ X, 
2 



where 



X = 



m 



A 2 \m 2 + M 2 



2 



By the help of the shift Y = X + X where 

v 1 (-2m 2 
X °~8\ 



3M 2 
2m 2 + M 2 



we obtain s = \Y ■ B ■ Y + s r 



The rotation matrix 



cos G sin G 
— sin G cos G 



(14) 



(15) 



(16) 



(17) 



with G = 7r/4 diagonalizes the quadratic form, 
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* = --^Y 1 ~ Y 2 ) 2 + \<X X + Y 2 ) 2 + s m . (18) 

The minimum is the result of the competition between the negative and the positive eigen- 
value, i.e. the trend to increase \Y l — Y 2 \ and to decrease \Y X + Y 2 \. The result is that 
the minimum is reached at the boundary of the quadrant X 1 > 0, X 2 > 0, i.e. there is no 
ferrimagnetic phase realized. Thus the mean-field vacuum is found at the minimum of the 
following two functions 

s(^ H , 0) = \m 2 <Pl + \<fa, s(0, ft H ) = l -C 2 {m 2 + M 2 )<f> 2 IH + \c^ IH . (19) 
The vacuum is antiferromagnetic when 



S 1 + M 

Ca \ m 



2 



X + — >! (2°) 



with the action density 



C 2 ( M 2 
] -CA 1+ m 2 



1 + 3 • (21) 



= 3 : The minimization of the action with respect 9 yields the condition (pH(f) IH sin 39 = 
0, whose solutions, <$h = 0, §m = and 9 = 0,7r/3 correspond to the anti-, ferro- and 
ferrimagnetic phases, respectively. The transformation 9 — > 9 + n/3, and <pH — » —<f>H leaves 
the mean field action invariant and it is sufficient to consider the case 9 = to explore the 
ferrimagnetic phase. The action density in the ferro- and antiferromagnetic phases is given 
by fllTf ) for 9^0. The antiferromagnetic phase is preferred for (|20| ) and the corresponding 
action density is (pT|) . For 9 = we have to minimize (|TTD numerically. 

N = 4 : The minimization with respect 9 yields (pj H sm49 = 0, showing the possibility 
of the ferrimagnetic phase when 9 = rwr/4. Since the action is an even function of the 
amplitudes 4>h and 4>ih for even iV we have the expression ATI) where 



*= l (nt f A „A B= X -( 1 (22) 

2 \C 2 (m 2 + M 2 )J ' 2 \3C 2 I J K ' 

what results 

1 ( m\l-3C 2 )-3C 2 M 2 \ 

X(9C 2 - C 4 ) V m 2 C 2 (C 4 - 3C 2 ) + C 2 C 4 M 2 J ' 1 6) 

The rotation (|T7|) satisfying the condition cot 20 = (1 — Ci)/QC 2 transforms the action into 
the diagonal form with the eigenvalues: 

^ (l + C 4 ± J(1 + C 4 ) 2 + 36C$-4C^ , (24) 

indicating that one of the normal modes is again unstable. Due to the negative eigenvalue 
and < < 7r/4, the minimum is always reached at the boundary of the quadrant X 1 > 0, 
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X 2 > 0, i.e. there is no ferrimagnetic phase. The condition for the antiferromagnetic phase 
and the expression of the action density are given by fl20|) and 

N > 4 : The procedure is the same as for N = 4, with the only difference is that C3 = 
and there is no ^-dependence in the sums C 2 and C 4 . 

Phase structure: The resulting action densities are summarized in Table |. The C4 de- 
pendence of the period length of the vacuum for c 2 = 2 and iV < N max = 32 is shown in 
Fig. [TJ. The general trend is that the increase of C4 pushes the minimum of the disper- 
sion relation p m i n towards zero thereby increasing the period length of the vacuum. At a 
certain threshold value the minimum at p m i n becomes so shallow that the potential energy 
turns the vacuum homogeneous and the system undergoes a ferromagnetic phase transition. 
Notice the usual signatures of the commensurate-incommensurate transitions, the "devil's 
staircase" structure. This is a competition between two length scales, the cutoff a and the 
period length of the condensate A mc , 

MA mc = Na. (25) 

The period length in lattice spacing units, iV(c 2 , c^)/M(c2, C4) is shown in Fig. [3] as the 
function of C4. It "locks-in", i.e. stays constant in a larger commensurate interval where 
the relative primes N and M are small. The long strips corresponding to the N = 2 and 
4 phases show a strong "lock-in" effect, contrasted with the gradual change of the period 
length for other values of N. The numerator and the denominator as the functions of C4 are 
non-monotonic in the same time. 

The phase structure in the plane (c 2 ,c 4 ) is depicted in Fig. 0. We identified the mean- 
field parameters N < N max = 32, M and cIaf for each point in the plane (c 2 ,c 4 ). The 
points where we enter into a phase with M = 1 by increasing C4 are indicated by the solid 
lines, the other phase boundaries are shown by dotted lines. The increase of N max makes 
the dotted lines denser without changing the solid lines or populating the white area. The 
vacuum in the upper left region is ferromagnetic. The lower right part of the c 2 , C4 > 
quadrant contains the inhomogeneous vacuua, each of them situated in slightly tilted strips 
with increasing N/M as we move upwards. The narrow triangular phase between the N = 2 
relativistic and non-relativistic phase is a relativistic phase with N = 3, see Fig. |^. It is 
interesting that the relativistic vacuua are realized for N = 2 and 3 only. We found no 
ferrimagnetic phase for the parameters considered. 



IV. EXCITATION BANDS 

The study of the more detailed structure of the dynamics starts with the determination of 
the possible excitations. The distinguishing feature of the antiferromagnetic vacuua is their 
inhomogeneity and the nonconservation of the momentum of the excitations. We introduce 
in this Section a convenient formalism for the description of the elementary excitations. 
Following the solid state analogies we rewrite our action by means of the Bloch waves which 
take the non-conservation of the momentum into account. 

Brillouin zones : First we split the first Brillouin zone 

B = {Pa.;M<^ A* = 1,2} (26) 
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of the d dimensional N antiferromagnetic phase into N dAF sub- zones within which the Bloch 
momentum is conserved. In the relativistic case when d^ F = d, we find 

K = {v»Av»-Pi N \*)\<^}i (27) 

where the center of the sub-zone, P( N \a), is given by 

^» = (28) 
in terms of the integer valued vector n^(a ) — 0, • • • , N — 1, 

a = l + J2M a ) N ^- (29) 

In other words, the integer component vector n^a) gives the center of the sub-zone B r a in 
units of 2tt/N. In the same time, it can be considered as an iV-base number. In this case 
its value labels the corresponding sub-zones. In the nonrelativistic case, dAF < d, we have 

K r = {p„ W - P { , N \ a )\ ^f' ^ d ^ 1^1 < > d AF ) ■ (30) 

In the next step we introduce a field variable in real space which is responsible for the 
fluctuations in each sub- zone, 

*><*> = L wr e "" m = L (S eW+P '" >W) < 31) 

by the help of the Fourier transform 

0b) = ^Ee-^ P 0W, (32) 

and its restriction into the sub-zones, <f> a (p) = <p(p + P^ N \a)). The computation what 
follows is considerably simplified if the N dAF Fourier transforms are extended over the first 
Brillouin zone as periodic functions, 

MP) = Mp+p {N) m, (33) 

where (3 is an arbitrary sub-zone index. The tilde on a momentum variable will always 
denote that the given momentum is in the sub- zone B\. The path integral is then written 

as 

n / dMOe-™ = nn / Mp)e- s[M , m 

p J a p J 

with 

1 r d d p 



LdS[<Pa] = 2 / ?2^ 0( ^ )[P2/C( " P2) + m2]m 



(2tt) c 

(2k)' 



+ I (A / f^^) ( 2 -)^(E^) (35) 
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where 



P»=4£sin 2 ^ (36) 

denotes the momentum square on the lattice. 

Flavor algebra : The N dAF sub-zones introduced above correspond to the different exci- 
tation bands of the coarse grained lattice whose lattice sites represent the primitive cells of 
the original lattice. The umklapp processes where a non- vanishing momentum is exchanged 
with the vacuum take the momentum from one sub-zone to another. To separate the sub- 
zone preserving and changing processes from each other we rewrite the momentum integrals 
in ( |35| ) as a sum over the sub-zones and integration within B%, 

r N d AF 

Jd d pf(p)= J2 J d d pf(p + P^(a)). (37) 

a=l 

The summation can be organized in a more transparent manner by considering the group 

Zd AF ,N = ® ff Z n (38) 
i=i 

which describes the shift of the momentum in the periodic, <1af dimensional region of the 
Brillouin zone. This group makes its appearance by considering the action (|35|) as a matrix 
element of an operator in the function space span by the "wave functions" <p a (p). It is 
advantageous to use the plane wave basis \P( N >(a) +p) = \a,p) in which the matrix element 
of the field operator 4> a (ft) is defined by 

< P', q\UP)\P, Q> =5 (pWffl + q + pW(a) +p- pW(?) - q') UP) 
= 5 (P W ifi) + PW (a) - (/?')) S (q + p- q') <P a (p) 

= (i a )^^Hq + P~q)MP), (39) 

where the periodicity (^) was used in the second line. The symbol 4> a (p) stands for operator 
when sandwiched between the basis vectors and for function in the c-number expressions. 
We introduced here a representation of Zd AF ,N by means of N dAF x N dAF matrices, 

d AF 

(l a ) P ,<r = II S <r lt +a ll ,p lt (modN) 1 ( 40 ) 
(t=l 

constructed in such a manner that 7" describes the effect of the exchange of the momentum 
p( N '(a) on the index labeling the sub-Brillouin zones. We shall need later a relation arising 
from the Abelian nature of the group Zd AF ,N, 

(r) a +p(modN),0+p(modN) = (Y^ (Y)^ 1 ) a,/3 = (^)a,/3, (41) 

where T is an arbitrary product of the 7 matrices. The action can now be rewritten in a 
more compact notation as 
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L d Si<P] = (O,O|i0(O)[P 2 /C(-P 2 )+m 2 ]0(O) + ^ 4 (O)|O,O) (42) 



where the operator 0(0) is given by 

*< > = /r£^> = ?/7Sd*<*>- (43) 



Elementary rearrangements yield 



L d S[<P] = tr{± J ^ d A-P)i>C(p)+m 2 m + \ (jlj |^tfP*)) W^E?*)}. (44) 



where $>= Q 7 a , and the matrices JC(p) and m 2 are given by 

£^(p) = * a /(P (ff) (a)+pm-? 2 (P^(a)+p)) 



m 2 aP = 5 a ^m 2 . (45) 



The action (44) includes iV dAF fields whose flavor mixing is handled by the matrices j a . In 
fact, the real space expression is 

S[<f>] = J d d xtr^-d^x)K(^^ny^x) + Vbtfx))}. (46) 

In lattice regularization the lattice sites correspond to the primitive unit cells of the original 
vacuum. Our goal, to trade the momentum non-conservation on an inhomogeneous vacuum 
into a multiplicity of excitation bands over a homogeneous vacuum, is completed. 

V. ELEMENTARY EXCITATIONS 

We determine in this section the elementary excitations which are the eigenfunctions of 
the second functional derivative of the action, evaluated at the tree level vacuum. We choose 
our tree level vacuum in the phase iV to be 

0(x)=^e- p(JV) M$ Q . (47) 

a 

We need the eigenvectors of 

G^(P) = f^f. r , =(%~)+m 2 + 3Af) (48) 

0<t>a{-p)0<t>p\P) |0=$ V Ja ^ 

The propagator can formally be written as 

N d AF 

G{P) = £ MP)K\P)£(P), (49) 

Q = l 

where 
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(£(p) + m 2 + 3Af ) Vo(p) = K(P)MP)- (50) 

The diagonalization of the quadratic part of the action for a given Bloch wave number 
provides the propagator and the corresponding band structure, A a (p). 

We start with small values of C4, i.e. we are either in the ferromagnetic or in the N = 
cIaf = 2 antiferromagnetic phase, when $ a = S a: i4>iH, or $ a = S a ^(j)jH, respectively. Since 
7 2 = 1 the inverse propagator is 



where 



and 



G- 1 (p) = m 2 MF +p 2 lC(-p 2 ), (51) 



p M = 2sin^ (52) 



2 _ 2 , ox j,2 _ / -2m 2 ferromagnetic, , . 

m MF -m +^i-|_ 2m 2_ 3 _ M 2 (1)2) ^ C2)C4) anti f erroma g ne tic. 

The propagator can be written as 

G- 1 = m 2 MF + V 2 (l- c 2 V 2 + c 4 P 4 ) , (54) 

or 

G a , f3 (p)=5 a , (3 G(P^\a)+p). (55) 
The long wavelength fluctuations give 



where 



and 



G~m = m 2 MF (a) + Z(a)p 2 + 0{p A ) (56) 



m 2 MF (a) = m 2 MF + M 2 (l, 2, m(a) + n 2 (az), c 2 , c 4 ) (57) 



Z(a) = l 1 , ,0 ! (58) 

v ; I -1 + 8c 2 - 48c 4 a = 4 v ; 

The inverse propagator always has a local minimum at "P 2 = 0. For certain values of c 2 
and C4 it may develop another minimum at V 2 = V 2 > when considered as a function of 
V 2 . This minimum is realized kinematically for V 2 < 8 only, when it is reached on a closed 
line in the Brillouin zone, c.f. Fig. |6|, a structure reminiscent of the roton spectrum ||. 
Thus the van der Waals-type force, represented by the choice c 2 > leads directly to the 
appearance of the additional minima of the dispersion relation interpreted as rotons [fLOfl . 

To follow this in detail, c.f. Fig. ^, we start with the condition for the extremum, 

G -1 (p) = 2 sinj9 M (l - 2c 2 P 2 + 3c 4 P 4 ) . (59) 



dp 
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Apart of the points p^ = Pj^'ia) which are always solutions, the root of the expression in 
the parentheses yields another extremum, 



so long as 



C4 < f • (61) 

The extremum at V 2 = V 2 is always a minimum. To select the other minima we need 
the second derivative matrix, 

y } = 25^ v cos Pll (l - 2c 2 V 2 + 3c 4 P 4 ) - 8 sm Pfi sinp„ (c 2 - 3c 4 V 2 ) . (62) 

The sub-zones a = 1 and 4 contain particle like excitations because the centers of these 
sub-zones are local minima, 

d 2 _ 

-G (p)| p=P (2) (1) = 2S^ U 



dp ii dp v 
2 

- — - — G^(p) \ P=P (2) (4 ) = -2^(1 - 16c 2 + 192c 4 ). (63) 

OPfiOPv 

In contrary, the other sub-zones contain saddle point only and do not support particle like 
excitations, 

d 2 d 2 

d 2 

= 'dp^ ' 1 ^^ 
d 2 

= -2(1 -8c 2 + 48c 4 ). (64) 

When c 4 is small enough then V 2 > 8 and the inverse propagator possesses two discrete 
minima, at p = p( 2 '(l) = (0, 0) and at p = p( 2 \4) = (it, ti). The absolute minimum is at 

_ J (0, 0) ferromagnetic, , > 

P \ (ir, 7r) antiferromagnetic. 

As the value of c 4 is gradually increased from zero V 2 reaches the value 8 and it is better to 
follow what happens in the two phases separately. In the ferromagnetic phase the degenerate 
local minima are found along a closed loop in the vicinity of the point p = (tt, it) when (60) 
reaches the allowed kinematical regime, V 2 < 8, for 

|<c 4 , and §-Y^<C4. (66) 
15 



This line becomes the absolute minimum when the smaller root of the expression 



1 - c 2 V 2 + c 4 P 4 (67) 
as the function of V 2 turns out to be smaller than 8, 

c 2 c 2 1 , , 

— < c 4 or C4 < -. (68) 

16 8 64 v ; 

In fact, flSTD is negative between the roots and the inverse propagator goes below its value 
at V 2 = according to (|54f) . 



It is shown in Fig. |3| that the line of degenerate local minima appears within the region 
bounded by the part ab of the parabola P, (|6lD and the line segments be and ac, defined by 
the first and the second inequalities of (|68| ) and (p6[) , respectively. The line of degenerate 
minima turns out to be absolute minimum below the parabole fl61~D and above the line 
segment be and the ferro- antiferromagnetic phase boundary. The absolute minimum of the 
inverse propagator is degenerate and lies on a closed line in either phase for c 2 > 1/4 when 
the inequality ( |6"BD is trivially satisfied. 

As the value of C4 is increased in the antiferromagnetic phase the inverse propagator 
keeps its minimum at p — (jr, ir) so long as the second inequality of (|66D is violated. The 
absolute minimum is degenerate and found along a closed loop around p = (ir, ir) when c 4 is 
further increased. We may avoid the phase N = 3 by choosing large enough c 2 and arrive at 
the N = 2 nonrelativistic phase without modifying the propagator at the phase transition. 
The value of the lattice momentum at the minimum, (0), is a monotonically decreasing 
function of C4 and reaches 4 for 

the line S of Fig. [|. At this point the degenerate minima of the inverse propagator form a 
square and for larger values of C4 it is deformed into a closed loop around p = 0. 

The further increase of C4 brings us to the higher N phase boundaries and the situation 
becomes more involved. 

We turn now to the question of the critical points. According to (^3|) and ( |57| ) the second 
order phase transition is reached when both m 2 — ► and M 2 — > for m 2 < 0. Let us take 

M 2 = -fi 2 M a 2 , m 2 = -fi 2 m a 2 . (70) 

The first equation and flTOD assert that the criticality is reached either for 

° 4 ^ 4cWsin 2 (f) ~ 164 F sin 4 (f ) ' 

or as N/M — ► 00. Both cases require the vicinity of the ferro- antiferromagnetic transition 
line, the value of c 2 and C4 is finite in the first case and diverging in the second. Since 
C±(N) > C 2 (N) (^l|) shows that [i 2 M > is needed to reach this phase transition. Thus the 
ferro-antiferromagnetic transition line is critical in the mean-field approximation. 

We close this Section with a remark about the Goldstone modes. The inhomogeneous 
vacuum breaks the external symmetries and we find Goldstone modes for the models in the 
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continuum. The lattice regulator reduces the space-time symmetries into a discrete group 
and there is no reason to expect massless phonons in the antiferromagnetic phases. But the 
strength of the breaking of the continuous part of the space-time symmetries is 0(M/N). 
Thus the continuous kinematical symmetries are restored as N/M — > oo and the Goldstone 
theorem is regained asymptotically, i.e. \{p m in) = 0(M/N). An important consequence of 
this gradual restoration of the Goldstone theorem is the absence of long range modes and 
the possibility of supporting the periodic vacuum in two dimensions when N is finite. This 
does not mean a long range order because the non-trivial part of the vacuum is squeezed 
within the distance Na/M — > 0. 



VI. ONE-LOOP RENORMALIZATION 

It is shown in this Section that the one-loop effective potential of our model can 
be made finite by the introduction of an appropriate mass counterterm and the ferro- 
antiferromagnetic transition remains a critical line on the one-loop level. 

It is straightforward to derive the Feynman rules for ([|4]) and to compute the radiative 
corrections. The one-loop effective potential, the generating function for the 1PI vertices, 



n=0 n - ai,-,a n 

can be written in the one-loop approximation as 

V eff ($) = V tree (®) + \J -^trHtip) + m 2 + 3Af ], 

where 

VtreeW = * (^(0)f + \m 2 f + ^ . 

We split the mass term into a renormalized and a counterterm, 

" - m 2 R + 6m 2 , 



(72) 



m 



(73) 



(74) 



(75) 



and make the replacement m 2 —>■ m? R in the loop integral. 

The ^-dependence of the one-loop integrals can significantly be simplified by using the 
relation (P) with T = 



(f)c 



: trf, 



(76) 



where no summation is made for the index a. Since the matrix K(p) + m? R is diagonal, 



V eff ($) = VtreeW 



d d p 

W) 



:tr In 



+ m R + 



3A 



tr(f) 



(77) 



The ultraviolet divergences are identified in two steps. First we expand in the field 
dependence and write the loop contribution as 
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i i #p r sxtrf y 

^2nJ (2 I r)" tr U(fi)+m|J ' { '*> 

to recover the usual one-loop resummation in the effective potential. In the second step we 
expand the integrals in the lattice spacing a, 



a 



p\<ir/N \p\<ir/Na 



X 



(79) 



We have at most a logarithmic divergence in two dimensions (n — 1) and the dimensionless 
inverse propagator, the integrand of the left hand side, has a finite, 1/a independent min- 
imum. Let us write the smallest eigenvalue of the inverse propagator around its minimum 

as 

Xa min (Prmn + P) = ™\ + Z p 2 + O(aY), (80) 

with Wq = a 2 m^j in , what allows us to identify the divergent part of the loop integral for 
d = 2, 



1 r d 2 f> 1 

' P (81) 



2 J (2nymi m + Z pi- 

\p\<-K/Na 

Thus the one loop ultraviolet divergence can be removed by setting 



|p|<7r/7Va 



w = - / (82) 



When = 0(a~ e ), e > 0, then no divergency arises. If there are several finite minima 

of the dimensionless inverse propagator then we sum over them in the counterterm. 

The continuum limit, ( [70]) can be achieved along the ferro-antiferromagnetic transition 
line of the plane (02, C4) . We have one or two particles by approaching this line from the ferro- 
or the antiferromagnetic phase below the point c = (1/4, 1/64), respectively. Above the point 
c the rotons appear in both phases and replace one of the particles of the antiferromagnetic 
vacuum. One can find specially interesting continuum limits in the vicinity of the point c. 
(i) Approach from below the line ab in the ferromagnetic phase, which gives a single particle 
and monotonically increasing inverse propagator with the momentum, (ii) The approach 
above the line ab but below be yields one particle and rotons. The roton momentum diverges 
with the cutoff, a reminescent of the new excitation bands in the antiferromagnetic phase. 
The value of the inverse propagator at the degenerate minima, the roton mass square is 
finite but larger than for the absolute minimum, at p = 0. Thus the rotons are heavier than 
the particle, (iii) Approaching above the line be yields a model with a particle and rotons 
where the rotons are lighter than the particle, (iv) The approach from the antiferromagnetic 
phase produces a model with two particles. 
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VII. NON-PERTURBATIVE VACUUM 



In order to assess the importance of the fluctuations, we performed a Monte Carlo sim- 
ulation of the model (f|). The resulting phase structure is summarized on Fig. [|. The 
coupling constants A = 0.05, m 2 = —0.1 were chosen and the simulation on the line C2 = 2 
was done at 200 x 200 lattice size. The ca dependence of the results was monitored with 
a particular care by scanning the region < C4 < 1 with step size AC4 = 0.01. The other 
points were obtained on 40 x 40 lattice. The performed tests have showed no appreciable 
finite size dependence in the qualitative and quantitative features of the phase diagram. We 
used the metropolis update algorithm, carefully checking the statistics to make sure that 
no statistical error would change our conclusions concerning the phase diagram. In order to 
test ergodicity, local cluster algorithms were tried as well as different initial conditions. 

The letters along the vertical lines of Fig. [| indicate the qualitative space-time struc- 
ture of the vacuum seen in the simulation: A=ordered (N = cIaf — 2) antiferromagnetic; 
L=labyrinths; W=plane waves; P=weakly antiferromagnetic (onset of the crossover on the 
finite lattice), F=ferromagnetic and F'=weakly ferromagnetic. To understand the phase 
structure we recall that the parameters of the bare action characterize the dynamics at the 
cutoff. Consequently one may extract useful informations about the short range order of the 
vacuum by considering the coefficients appearing in the quadratic part of The sign of the 
coefficients determines the ferro- or antiferromangetic nature of the couplings. We introduce 
the subscript (n, m) which shows the separation of the two field variables multiplied by the 
given parameter, v4 ,o = 2 + m 2 /2 — IOC2 + 56c4, -81,0 = — 1 + 8c 2 — 57c 4 , 6*2,0 = — c 2 + 12c4, 
Ei t i = — 2c2 + 24c4, -D3.0 = — C4 and i*i i2 = — 3C4. The lines where A, B, C and E change 
sign, 

10 m 2 2 . 4 , 

Ca = Co (A), 

56 112 56 v 1 

C4 = ^c 2 -- (5), 

c 4 = ^c 2 (C,E) (83) 

are shown in Fig. |] by solid lines. These lines intersect and the short range order varies 
in a complicated manner for c 2 < 0.45. But once this value is reached, the sequence of 
the change of the signs as C4 increases is always the same. The importance of these signs 
is that they create frustrations whenever C4 > 0. There are no frustration when ca = 
because D 3 = F 21 = and all others favor the N = 2, d^F = 2 antiferomagnetic or 
the ferromagnetic vacuum for c 2 > 1/8 or c 2 < 1/8, respectively. For c 4 > but below 
the line C, E, in the antiferromagnetic N = 2, (Iaf = 2 phase B 10 > 0, C 2)0 , £a,i < 
favor this kind of vaccum. But the other signs, D^o, ^2,1 < introduce frustrations whose 
density increases with C4. The further increase of C4 flips the sign of B and ultimately A, 
destabilizes the lock-in mechanism at iV = 2 and opens the way for the rapid variations of 
the devil's staircase. In between the lines C, E and B, C2 0, £1,1 > and the frustration 
density is increased because only 5 10 favors this vacuum. The result is a strong increase of 
the fluctuations. The mean-field approximation is obviously unreliable in this regime and 
the simulation produces labyrinth-like vacuum, see below. For large C4, in the ferromagnetic 
phase -Bio, D 30 , F 1 2 < act in favor of the homogeneous vacuum but C 2 o, En > 
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generate frustrations which might explain the weakening of the ferromagnetic condensate 
for large C4. For the intermediate values of C4 the competition between the different terms 
is more involved and the compromise between the different competing terms is reached over 
a longer length scale according to the mean-field solution. The qualitative conclusion is the 
separation of the stable low N = 2, cLaf = 2 vacuum from the high C4 ferromagnetic and 
an intermediate C4 strongly frustrated antiferromagnetic vacuua. 

The frustrations are the lattice defects of the antiferromagnetic vacuum. In order to 
understand their production mechanism in the N = 2 antiferromagnetic phase we used two 
different initial conditions, an ordered and a disordered one, 



Lord . 
init ' 



x)=M-ir +x , (84) 



and 



idis 
init 



x) = 0, (85) 



respectively. The initial condition 4>f n s it ^ x ) which actually looks ordered becomes disordered 
after few Monte-Carlo sweeps. This happens because it represents an unstable equilibrium 
position and the local field variable "rolls down" to one of the potential minima (j){x) = ±4>ih- 
This leads to the formation of iV = 2 antiferromagnetic domain structure. These domains 
are separated by walls of links where the field variable has the same sign. The region on 
the plane (02, C4) where the antiferromagnetic domains consisting of the patches <j)°^f t (x) and 
—(f)°nf t (x) develop long and winding boundaries is called labyrinth-like and is denoted by L 
in Fig. f§ 

For large values of C2 the amplitude of the modulated vacuum is large so the frustrations 
move very slowly in the simulation time. The antiferromagnetic domain walls turned out 
to be very slow variables, as well. The domain walls were always generated below the line 
C, E of Fig. f| when the disordered initial configuration was used. The vacuum obtained 
by the runs with ordered initial configuration did not support the domains. The question is 
whether the domain walls are real degrees of freedom or reflect the insufficient convergence 
of the simulation method. We have developed cluster algorithm and found that for small 
C2, (C2 < 1-0), the walls have dynamics, reach an equilibrium and may disappear. For larger 
C2 the domain wall motion slows down despite of the cluster algorithm. The thermalization 
was safely reached within each domain. It remains an intriguing question if the domain walls 
thermalize in this regime with extremely long relaxation time, i.e we are in a glassy regime 
22" I or the ergodicity is definitely lost and the vacuum consists of a stable, disordered network 



of domain walls. In this regime the frustrations act as scattering centers without feedback 
from the fast dynamics of the elementary excitations, a dynamical situation reminiscent of 
the quenched disorder in solids. For C2 < 1 the mean field value, 4>ih, is small enough to 
make the domain wall fluctuations more likely and the cluster update averages over the 
different rearrangements of the domains. 

The different regions shown in Fig. |] were studied in more detail at C2 = 2. The negative 
action density obtained by starting at C4 = and sweeping the interval < C4 < 1 is plotted 
in Fig. |5|(a). The results obtained by the ordered and disordered initial configuration are 
indicated by + and squares, respectively. The ordered vacuum has lower action density up 
to C4 ~ 0.2. The mean-field action density is shown by a solid line. It is instructive to 
follow the second lowest mean-field solution, indicated by the circles. The splitting between 
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the lowest and the second lowest mean field action level reflects the stability of the vacuum 
against the change of the long range order. The agreement between the Monte Carlo and the 
mean-field results is remarkable for small c 4 in the relativistic phase. Right at the relativistic- 
non relativistic phase transition, the ordered configuration is a rather poor approximation 
and the disordered vacuum adjusts itself easier to the value cIaf = 1- For c 4 > 0.3 the 
mean field badly over estimates the true vacuum action density, indicating the presence of 
strong fluctuations. Notice that the phases with N = 2 and 4 are more stable against the 
modification of the long range, in agreement with their stronger "lock- in" in Fig. [[]. 

Another insight into the vacuum can be gained from the inspection of the period length 
of the vacuum, A mc , measured by 

< 2 > WvacJ 

The numerical and the analytical results depicted in Fig. |5](b) show that the first order 
transitions of the mean field approximation between the relativistic and the non-relativistic 
N = 2 antiferromagnetic phase is smoothened out in the simulations when the disordered 
initial condition is used. The ordered initial condition follows the mean field curve within the 
phase N = d^p = 2. The slightly higher action of the disordered initial configuration runs 
indicates that the true vacuum is close to being ordered and the relativistic-non relativistic 
phase transition is of strongly first order. One is tempted to conclude that the fluctuations 
smoothen out the commensurate-incommensurate transitions but better statistics is needed 
to settle this question in a satisfactory manner for the whole phase diagram. The difference 
between the numerical results and the mean-field solution is the largest in the phase N = 
2, cIaf = 1- The mean-field approximation slightly underestimates the period length of the 
vacuum in the vicinity of the ferromagnetic transition. This is consistent with one of the 
remarks made about Fig. |](a), namely that the fluctuations in this regime are stronger 
than expected by the mean-field approximation. In fact, the stronger fluctuations lower the 
critical value of c 4 so the period length of the vacuum diverges faster in function of c 4 than 
in the mean-field expression. 

The strength of the modulation of the vacuum is displayed in Fig. |5](c). It is simplest to 
express it in terms of the Fourier transformed field 

m = ^<P{x)e- vp - x , (87) 

obtained on an L x L lattice. We may split the expectation value < |0(p)| 2 > into the sum 
of the condensate and the fluctuations, 

< |0(p)| 2 >=< |0(p)| 2 > c + < |0(p)| 2 > fl , (88) 



< \4>(p)\ 2 > c = I < }(p) > | 2 , < \4>(p)\ 2 > fl = G c (p). (89) 

were G c (p) stands for the connected propagator, given by (|54]) in the leading order of the 
perturbation expansion. A simple approximation for the strength of the modulation of the 
vacuum is 
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4> 2 C = max p (\(j)(p)\ 2 ) c ~ max p (\<fi(p)\ 2 ) ~ (max x (j)(x)) 2 ~ (min x (j)(x)) 2 . (90) 

max p (|0(p)| 2 ) is displayed by plus and square for the ordered and disordered initial config- 
uration, respectively. The star shows the value of the last two expressions in (^). We find 
that the different estimates for the strength of the condensate agree in the relativistic phase 
except the simulation results corresponding to the disordered start. The difference between 
the mean-field solution indicated with x and the numerical results can be considered as a 
measure of the strength of the fluctuations. Note the local maximum and minimum in the 
fluctuations at C4 ~ 0.35 and C4 ~ 0.6, respectively. The former is in agreement with the 
remark made for Fig. [|(b). The average of the extrema of the field in the real space is higher 
than the mean field value and shows no structure as the function of C4 indicating that the 
fluctuations tend to lower the local values of the field variable |0(x)|. As noted before, the 
fluctuations increase as the ferromagnetic phase is approached, a result consistent with the 
second order nature of the ferromagnetic phase transition. 

The numerical results for the complete propagator, (^) what displays the structure of 
the vacuum and the elementary excitations in the same time are presented in Fig. |6|. We 
divided the interval [0, max p (\(j)(p)\ 2 )} into five equal segments and their contourplots are 
shown in the figures. The strength of the contour line increases with the amplitude, so the 
blacker regions of the plots indicate the location of the maxima. 

< Cd < 0.16: When c 4 is small we are in the relativistic N = 2 phase and < |0Qo)| 2 > 
depends strongly on the initial condition of the Monte Carlo simulation. The ordered initial 
configuration fl84|) yields a single peak for < |0(p)| 2 > at p = P^ 2 \A) = (tt,tc) suggesting 
little disorder. In the case of the disordered initial configuration, ^fnitt one nn ds a domain 
structure and the Fourier transform | < <fi(p) > \ 2 shown in Fig. ||(a) at C4 = is depleted at 
p = p( 2 )(4), < |0(p( 2 )(4))| 2 >m 0, and develops a ring of maxima around this momentum. 
It is easy to understand the minimum at the center. In fact, assuming that the vacuum 
consists of the domains of 0°^(x) and —4>inu( x ) i n ^ ne fractions c and 1 — c of the volume, 
respectively one finds 2 on(J = 2 H (1 — 2c) 2 /4. The domain pattern develops after few sweeps 
and the domain walls turn out to be rather slow variables. 

0.16 < Cd < 0.22: The excitations become more involved in this regime. The coefficients 
C and E of the lattice action are positive for C4 > 1/6 making the frustration density 
higher. Furthermore the propagator develops a circle of degenerate maxima around the 
point p = P^ 2 \A) for C4 > 31/192 m 0.16. The result is a cusp in the condensate as the 
function of C4 cf. Fig. |5](c), and the softening of the modes giving an increase of < |</>(p)| 2 >. 

0.22 < Cd < 0.3: The ordered initial configuration simulation recovers the right vacuum 
at C4 ~ 0.23 in a discontinuous manner, c. f. Fig. |6|(b) and (c). For C4 > 0.23 the simulations 
corresponding to the two different initial conditions yield the same result. The hysteresis in 
the C4 dependence, i.e. the later appearance of the ring for the ordered initial configuration 
case compared with the unstable starts suggests that the transition cIaf '■ 1 < — ► 2 is of 
first order. The roton minimum in the dispersion relation tends to break the straight lines 
where the frustrations are found and to distribute them in a more spherically symmetrical 
manner, resulting in more disorder and creating a labyrinth structure instead of the ordered 
nonrelativistic antiferromagnetic vacuum. 

0.3 < Cd < 0.38: The minimum of (Q) is the longest, being a square, for c 4 = 0.31 
according to flBTf). A typical example shown in Fig. |6](d) witnesses that the fluctuations are 
the strongest in this regime, when the volume of the phase space occupied by the soft modes 
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is the largest. As the length of the roton minima starts to shrink for ca_ > 0.31 so does the 
strength of the fluctuations as seen in Fig. ||(c). 

0.38 < Ca < 0.6: The period length of the vacuum growths beyond 2 at c 4 ~ 0.38 and 
< |0(j»)| 2 > becomes strongly suppressed at p = (0, n) and p = (it, 0) according to Fig. |](e)- 
(f). The rapid variation in C4 supports the discontinuous nature of the transition N = 2 — > 3. 
The further increase of C4 makes the restoration of the rotational symmetry more difficult. 
This is because the symmetry restoration is achieved by forming domains where the plane 
wave modulation of the vacuum has different orientation. The longer period length of 
the vacuum makes more energy consuming to break the straight plane wave by changing its 
direction, i.e. the domain wall energy density increases. This explains the breaking up of the 
ring into smaller segments and its ultimate reduction to few discrete peaks, as demonstrated 
in Fig. |(g). 

ca ~ 0.6: The fluctuations reach a minimum in the middle of the "lock-in" interval 
N/M = 4 seen in Fig. |l| The point where N/M is approximately the power of the 
highest order terms in the field variable of the lagrangian the above mentioned decrease 
of the fluctuations comes to a halt and is turned into an opposite trend due to the approach 
of the ferromagnetic phase transition. In fact, as the length of the roton minima shrinks 
and becomes less important the increase of C4 ultimately leads to the disappearance of the 
condensate which is triggered by the increase of the fluctuations. This is due to a tree-level 
effect, C4 makes the modulation of the vacuum more energetic so the amplitude of the mod- 
ulation decreases with the increase of C4. As the amplitude decreases it becomes easier to 
break a plan wave into domains what amounts to the attempt for the restoration of the 
rotational symmetry. The maximum in < <p{p) 2 > spreads from a well localized point over 
the whole circle of the rotons as C4 is increased beyond 0.6. 

0.6 < Ca < 0.98: As C4 increases and p m in approaches zero the dominant fluctuations are 
grouped on a circle around p = with increasing strength as shown in Fig. |6|(h). The 
condensate weakens and increasing period length in lattice spacing units is in agreement 
with the one-loop renormalizability established in the previous section. 

ca ~ 0.98: The precursor of the transition to the ferromagnetic phase is the appearance 
of a peak in < |0(j»)| 2 > at p = for C4 ~ 0.98. The further increase of C4 brings us into the 
ferromagnetic phase with roton excitations. 

We close this section with a remark concerning the Lifshitz point. It is a tricritical point 
at g% = c\ = in the model fl5|) where the wave vector of the periodic vacuum tends to 
zero as C\ — > 0~. This is to be compared with the "Lifshitz line" of our model, the curve 
separating the ferromagnetic and the antiferromagnetic phases in Fig. |2|. For a given C4 > 
this line gives a critical point which is reached by tuning C2 and m 2 . Since the coefficient of 
the lowest order term of the gradient in the action is kept constant [c\ = 1) the fine tuning 
of the higher order coefficients generates discontinuity for the wave vector of the vacuum at 
the critical pointQ. The peculiarity of the extension of the Lifshitz point to a line is that 
the wave vector of the vacuum is a discontinuous function of the coupling constants either 



1 The analogous situation at the usual tricritical point for the model (||) with C2 = 0, g§ > is the 
first order phase transition in the function of 34. We do not need C4 > in our case because the 
anharmonic term of the lagrangian stabilizes the vacuum for C2 < 0, C4 = 0. 
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when one crosses the line or when one moves along it in the antiferromagnetic side. One 
point of the Lifshitz line shows a further interesting feature, the even and the odd sublattices 
decouple at = l/4d, C4 = 0. This can be understood by checking the invariance of the 
propagator under the replacement V 2 —> 4d — V 2 , and allows us to construct the continuum 
limit of chiral bosons ||. 



VIII. WICK ROTATION TO REAL TIME 



Our model (P was obtained by the elimination of some particles and the truncation 
of the resulting effective action. The usual strategy to avoid the instabilities in real time 
mentioned above is based on the introduction of new, auxiliary fields which render the kinetic 
term quadratic. The most natural attempt is suggested by the formal similarity between 
the Pauli-Villars regulated models and theories with higher order derivatives fll3f , 0. In 
fact, by assuming that the inverse propagator (^) as a function of p 2 has simple roots only, 
we can write the propagator as 

(*(P ) ~~ /„? 1 !\/r2\/„9 1 !\^2\/„9 1 j\/t2\ ~~ „9 1 T72 „9 i TFi Z3 i T72 ' (91) 



where 



(p 2 + Mf)(> 2 + M|)(p 2 + Mf) p 2 + Mf p 2 + M$ p 2 + Mf 



Z- 1 = (M 2 2 -M 2 )(M 3 2 -M 2 ), 

Z- 1 = (M 2 - M 2 2 )(M 3 2 - M 2 2 ), (92) 



Z' 1 = {Ml - Ml){Ml - Ml 



.3 . 



(93) 



Since 



Z^Z^Z^ = -(Ml - Ml) 2 (Mi - Mf) 2 (Ml - Ml) 2 , (94) 



one of the coefficients Zj is negative when the roots M 2 are real, i.e. in the case of the 
antiferromagnetic vacuum. We have two complex roots, say M| = M|* in the ferro- or the 
paramagnetic phase. This gives Z 2 U = Z^ 1 and 

Z 2 l + Z, 1 * = (Ml - Ml) 2 < 0, (95) 

and makes ReZ^ = ReZ^ < 0. Thus the real part of at least one of the contributions in ( pl|) 
is always negative. 

Such a propagator can formally be obtained as the functional derivative of a Gaussian 
generator functional, 

3 

Z [j(x)} = / Dtyfe Ci , (96) 



where 

G jik (x,y) = -5(x - y)5 hk Zj 1 (U y - M 2 ). (97) 



24 



Note that the convergence of the functional integration requires that the field 



with 



ReZj < be purely imaginary, an indication of the presence of the negative norm states 
in the canonical quantization procedure. The partition function of the model ([1]) can be 
written as 



D[<j>]e 



■J d d xV{ 



s 



Zoij] 



(98) 



in the framework of the perturbation expansion. The Wick rotation of this path integral is 
straightforward and leads to a model with three particles and negative norm Hilbert space. 
The positive and the negative norm particles are mixed by the vertices of the interaction 
V(<f)). In order to have a physically acceptable model we have to ensure that the time 
evolution remains unitary when restricted to the subspace of positive norm physical states 
and the Hamiltonian is bounded from below. 

The unitarity of the S-matrix can be established for energies below the threshold of the 
negative norm particle production, the key element for the applicability of the Pauli-Villars 
regularization. The stability of the vacuum is reached in this regularization scheme by 
first performing the renormalization in the Euclidean space and making the Wick rotation 
back to real time after that. The negative norm regulator particles are suppressed during the 
removal of the cutoff thereby the Wick rotated renormalized theory possesses a stable ground 
state. The non-commutativity of the renormalization and the Wick rotation can simply be 
understood by noting that the formally Lorentz invariant Pauli-Villars regularization scheme 
violates a simple rule: The non-compact nature of the Lorentz group excludes the Lorentz 
invariant regulator schemes. 

What we have shown so far is that our model yields an acceptable theory for real time only 
after the perturbative renormalization in Euclidean space-time. One wonders if canonical 
quantization which leads to a non-perturbative formalism can carry us further in establishing 
an acceptable, non-perturbative theory. In our case it amounts to the introduction of the 



generalized coordinates 



n 



rii 



5L 
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6dM~ 



0, 1, 2, and momenta 
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SL 




(99) 



One can easily check that the Hamilton- Jacobi equations of the Hamiltonian 



H 



I 



d 



d-l. 



X 



EM 



Wj 



dd(J), ddd(f>) 



(100) 
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-hti + 0oV 2 0o) - c 2 ^r(0 2 - V 2 o ) 2 



2 v " i ' " u ^ z 2A 2 

- V 2 o )V 2 (0 2 - V 2 o ) + ^(0o) 



c 4 /(2vr) 2 \ 2 



2 V A 2 J 

are equivalent with the Euler-Lagrange equation. 

The quantization procedure is based on the canonical commutation relations 

6(x°-y°)[^(x),U k (y)} = i5 jtk 5(x - y). (101) 

The inner consistency requires that X and Hi be antihermitean operators. The standard 
path integral representation [E3| results 



^0,-01,02 



e" iTH |0o, 0i, <h) = J L>[0] J D[n]e^o Td *{E J n ^o^-^[n^]} ) (1Q2) 



where the fields Hi(x) and 0i(x), as the eigenvalues of antihermitean operators are purely 
imaginary. This feature, the negative norm nature of the particle corresponding to the 
auxiliary field 0i = <9o0 makes the real time phase space path integral ( |102j ) divergent and 
ill defined because the Hamiltonian (|100|) is complex, 



h[Uq, iii, n 2 , , </> u 2 ]* = h[Uq, -iii, n 2 , 0o, -0i, 02] 

^^0,^,^,00,01,02]. (103) 

In view of these failures of rescuing the model by the introduction of the auxiliary fields it 
seems remarkable that the Wick rotation based on our lattice regularized Euclidean model 
does yield an acceptable theory when the elementary fields are chosen with more care. 
We shall show that the cutoff theory is unitary and stable in the Minkowski space-time. 
The Lorentz symmetry is recovered in the renormalization process only, when the cutoff is 
removed. 

We start by noting that the elimination of a heavy particle which gave rise the higher 
order derivative terms in the effective action actually creates mixed states from the pure 
ones. Thus the effective theory should be recasted in terms of its density matrix. This is 
not what happens in the usual blocking procedure which yields the path integral expression 
as a representation of the matrix elements of the time evolution operator between pure field 
eigenstates. But the initial and final field configurations of the path integral for the effective 
theory specify the states of the light particles only. The information concerning the states 
of the heavy degrees of freedom is lost in the elimination process, by the trace operation in 
their Hilbert space. In other words, the initial and final states of the effective model specified 
in its path integral are mixed. The additional correlations in the expectation values arising 
from the off-diagonal elements of the density matrix are represented by the higher order 
derivative terms of the effective action. 

Suppose that we truncate the effective action up to 0(d 2n ). The correlations, i.e. the 
information about the heavy particles in the effective theory extend up to n lattice spacing. 
The states appearing in the path integrals can formally be considered as pure ones if these 



26 



correlations are properly kept in them. This can be achieved by the introduction of n — 1 
auxiliary fields at each site, like with the Pauli-Villars fields (|96| ) or the canonical quantiza- 
tion ( |102| ). Another possibility is the regroupment of n consecutive lattice sites in the time 
direction with their field variables into a single rz-component fieldF], in a manner similar to 



the one employed in ref. [25]. We follow the latter strategy, the reminiscent of the spin-flavor 



assignment for the staggered fermions in real-space time ||26|| . The corresponding action for 

4 j (t,£) = <i>(7a + j,x), (104) 

where t is integer and j — 1, • • • , n is 



(105) 



t,x,y 



t.x 



where the space-time vectors x = (t, x) label the 'fat' lattice sites with n degrees of freedom. 

In the case of the action (|J) we regroup three time slices into one new 'fat' time slice 
(n = 3), (f>j(t, x) = 4>(3t + j, x) for j = 1, 2, 3 and the matrix A"^ connects the (j, k) pairs 
shown in Fig. [7], 

A~l(x, y) = 2B5 Si g5 jt i5 k) 3 + 2Cfe^(5 ii2 4,3 + ^,1^,2) + 2D5g^S jtk 



(106) 



in two dimensions (F = 0). The transfer matrix for the field <f>j is given by 

T = e~ 3aH (107) 

in terms of the original Hamiltonian. 

Our reasoning is the following. The unitarity of the time evolution operator for the states 
created by the field operator <p(x) will be the result of the realness of the spectrum of the 



square of the transfer matrix with n = 3, T 



-6aH 



The positivity of the spectrum of T 



and the positive definiteness of the norm in the Hilbert space will be proven by verifying 
that the reflection positivity holds for T [2"3|. The cancellation of the runaway amplitudes 
constructed perturbatively in ref. [14| is assured by the reflection positivity for T in a non- 
perturbative manner. The stability of the vacuum, i.e. the boundedness of the spectrum 
of the Hamiltonian from below follows from the finiteness of the regulated Euclidean path 
integral. 

The proof of the reflection positivity relies on the Euclidean time inversion operator 
defined as 



2 The need of the increase of the number of degrees of freedom can be understood in the semi- 
classical limit, as well. The unique characterization of the trajectories interpolating between the 
initial and the final field configuration requires the knowledge of 2n quantities. 
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Q(t,x) = (-*,£), e[(f)j(x)} = 4 -i(Oa;), (108) 

where the transformation j 1 — > 4 — j performs the time inversion for the degrees of freedom 
within a 'fat' lattice site. Note that the base point of the inversion on the original lattice 
for a given n is at a lattice site or halfway between two such sites for n odd or even, 
respectively. The time inversion B acts in the external, space-time coordinates and the 
internal space of j or other indices, it changes the sign of the time coordinate and performs 
the necessary inversions in the internal space, such as the transformation j — > 4 — j, the 
complex conjugation in the case of complex quantities, respectively. It is useful to introduce 
the internal time parity Tp — ±1 as the eigenvalue of the time inversion in the internal space 
of an observable F satisfying the condition 

[F[(f)(x)}} = T F F [<f>(Bx)] . (109) 

The reflection positivity of the transfer matrix T can conveniently be expressed by means 
of the time inversion O, it is the requirement that the inequality 

(F6[F])>0 (110) 

hold for any local functional F[<f>] of the field variable taken for t > 0. The time evolution 
is clearly unitary in the subspace of the Fock space which is span by the states obtained by 
the application of the operator F[<p} on the vacuum. 
To prove ( |110|) we write the action ( |105|) as 



S = S + S- + S + (111) 



where 



s = E^(<MM)), 



S-= E Mt + l,2)^ k {x,y)<t>k{t,y) + £ V&j^x)), (112) 

t<0,x,y t<0,x 

s+= E ^(t + M)A^(f,y)</> fc (t,f) + e y(4>j(t,x)). 

t>0,x,y t>0,x 

We shall need below the relation 

e[s ± [<p]] = s T [<p], (us) 

a consequence of the time reversal invariance. In terms of the matrix Ajl(x, y) this amounts 
to the equations 

e[ E J (t + i,x)ATi(f,y)0 fc (t,f)]= ]r e[^(t + i,f)]e[A^(f,y)]e[</» fc (t,x)] 

t>0,x,y t>0,x,y 

= E 04-i(-t-i,f)e[A-i(f,y)]</) 4 _ fc (-t,f)] 

t>0,x,y 

= E 04- fc (t + l,x)e[A^(x,y)]0 4 -,(t,x)] 
t<0,x,y 

= E Mt+l,x)e[A 4 l jA _ k (x,y)}^(t,x)} 

t<0,x,y 

= E <^ + l,x)A^£, $</>,•(*,£) (114) 

t<0,x,y 
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which yield immediately 

Q[Aj t l(x , y)} = Ai-Vite = A^(f , f). (115) 
The left hand side of the inequality (|110| ) can now be written as 

(FQF) = [ D t=0 [4>]e- S °^ J A<o[0]e- s - [ ^ )] 6[F[0(x)]] J D t>0 [<p}e- s ^F[^x)} 
D[(f)]e- SoWx)] Q [ e - 5 +[^)l] Q[F[(J)(x)]]e- s + Wx)] F[(J)(x)] 
D[(f)}e- SoWx)] e- s+w@x)] T F F [0(0x)] e - s+[,K * )] F[(p(x)], (116) 



where we assumed in the last equation that the functional F possesses a definite internal 
time inversion parity, Tp = ±1. This allows us to eliminate from the integrand by writing 

(FQF) = J A=o[0]e- 5o T F (J A>o[0]e~ 5+ F[0(x)]) 2 . (117) 
The result is the inequality 

T F (FQ[F]) > 0, (118) 

holding for any real functional F. 

Let us denote by Hi the subspace of the Fock space 7i which is span by the application 
of the +1 internal time parity local functionals of the field operator 4>(x) on the vacuum and 
by H rth its orthogonal complement. We have 

H = H 1 ®H orth (119) 

in an obvious manner. The time evolution is unitary within TCi, and the vectors in 7i or th, 
among others the ghost states of ref. [@, have negative norm in agreement with ( |103| ). 

More care is needed in the presence of a time dependent condensate. This is because 
the condensed mode is treated classically and the problems of the quantum treatement, 
in particular the unitarity of the time evolution and the reflection positivity apply for the 
quantum fluctuations only. The reflection positivity should hold for the dynamics of the 
quantum variable <j> q (x) = 4>(x) — (4>(x)) which is governed by the action 

SM = S[cf> q + (cf>)]. (120) 

Note that ( |1 13| ) is satisfied by S q , 

Q[s±[<p q + (m = s T [<p q + m, (i2i) 



according to ( [L15|) . This relation allows us to repeat the steps ( |116| )-( pT7|) after having 
performed the replacement S[<f>] — > S q [<f>]. 

Another problem one has to take care is the choice of n, the number of degrees of freedom 
at a 'fat' lattice site, n must be at least half of the highest power the gradient operator 
appears in the action in order to decouple the next-to-nearest neighbors. In the same time 
it must be integer times the period length of the elementary cell in the antiferromagnetic 
phase. These conditions require the choice n = 4 when N = 2 and n = N for N > 3. 

The reflection positivity satisfied in the subspace Hi of the N = 2 antiferromagnetic 
phase indicates that one has to use the space parity eigenstates, the superpositions of the 
two chiral bosons observed at c 2 = l/3d, c 4 = |, because only the scalar particle is 
physical, the pseudoscalar modes correspond to ghost states. 
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IX. CONCLUSIONS 



A simple two dimensional effective theory was considered in this work for an elementary 
interaction which is strongly repulsive for short distances and attractive at long distances 
and supposed to form a crystalline ground state. Such an interaction is coded in terms of 
the effective action with higher order of the derivatives. The inhomogeneous ground state 
is reproduced in the saddle point approximation of the effective theory. 

The mean-field, tree level phase structure of the model is quite involved and displays 
several inhomogeneous phases with a number of commensurate-incommensurate transitions. 
The elementary excitations contain modes similar to the rotons of liquid He 4 but the exci- 
tation spectrum is usually massive in lacking of a conserved particle number. Our effective 
theory is a simplified version of the higher dimensional systems with inhomogeneous ground 
states, such as the solids, antiferromagnets and materials with charge density phase. The 
formation of the inhomogeneous vacuum by the condensation of particles with non-vanishing 
momentum is a manifestation of the additional soft modes which characterizes the quantum 
phase transtitions. 

The one-loop corrections give a line of ultraviolet fixed points with variable particle 
content. The period length of the vacuum can thus be sent to zero. The Poincare symmetry is 
restored in the continuum limit and the vacuum becomes homogeneous for the measurements 
made at finite energies. In the same time the excitation spectrum of the model remains 
nontrivial, reflecting the inhomogeneity of the vacuum. The different dispersion relation 
branches are interpreted in the continuum limit as excitations with different flavor. This 
unusual vacuum-excitations correspondence opens the way for the construction of new kind 
of unified quantum field theoretical models where several particles are described by the 
same field. Our effective theory was found to have a consistent extension to real time in the 
even internal time parity sector of the Fock space despite the presence of the higher order 
derivative terms in the action. 

The numerical analysis performed by the Monte Carlo method confirms the mean-field 
prediction of the phase structure and is consistent with the criticality at the ferromagnetic 
phase transition line. The frustrations are the slow modes of the simulation when the 
amplitude of the modulated vacuum is large suggesting the possibility of a quenched disorder 
variable in real time, as well. 

We note finally that there are no massless excitation modes above the modulated vacuum 
with finite period length in lattice spacing units. This is specially striking in two dimensions 
where the periodic structure of the periodic vacuum is not necessarily destroyed by the 
Mermin-Wagner-Coleman theorem [p7f| . The continuum limit of the antiferromagnetic phase 
where N/M — > oo might be similar to the planar X-Y model with power like decay of the 
correlations. 
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TABLES 



N = 2 
N = 3 
N = 4 
AT > 4 



qUf = l 

T-^ I 4(l-4c 2 +16c 4 ) 

2 /, , 3(l-3c2+9c 4 ) V 

3 V 1 + ^ J 

(1 + 2(1 -2c 2 ) 2 

2 /i , 4s 2 (l-4s 2 c 2 +16s 4 c 4 ) V 

3 y J 



c/af = 2 



'j , 8(l-8c 2 +64c 4 ) 



1 + 



6(l-6c 2 +36c 4 ) 



+ 4(l-4c 2 +16c 4 ) \ 



1(1 + 



2 (l-8s 2 C2+64s 4 c 4 ) 



TABLE I. The antiferromagnetic action density, s(c/af, AT)/s(0, AT). The phase angle has been 
chosen to be a = tt/4 for N = 4 to minimize the action and s = sin7r/A^, N > 4. Since s(0, A") < 
the antiferromagnetic phase is preferred against the ferromagnetic one whenever the corresponding 
expression in the table is larger than 1. 
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FIGURES 

J I L 
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C4 

FIG. 1. The devil's staircase, the C4 dependence of the period length of the vacuum at C2 = 2. 
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0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



C2 

FIG. 2. The mean field phase structure in the (02,04) plane. F denotes the ferromagnetic 
phase. The lower edges of the phases with M = 1 and M > 1 are indicated by solid and dashed 
lines, respectively. The values of N, cIaf are shown for M = 1 at the right. 
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0123456789 10 

C2 

FIG. 4. The phase structure and the frustrations on the plane (02,04). The lines where the 
coefficients A, B, C and E of the lattice action change sign in the (02,04) plane. The letters 
along the vertical lines are to indicate the qualitative space-time structure of the vacuum seen in 
the simulation: A=ordered N = 2 antiferromagnetic; L=labyrinths; W=plane waves; P=weakly 
antiferromagnetic (onset of the crossover on the finite lattice), F=ferromagnetic and F'=weakly 
ferromagnetic. 
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FIG. 5. Averages at C2 = 2, (a): the action density, — < S > /L 2 , (b): the wavelength of 
the modulation of the vacuum, measured by — < > / < 4> 2 >, (c): the amplitude of the 

modulation. 
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(d) (e) (f) 




FIG. 6. The Fourier transform < |0(p)| 2 > in the plane (pi,P2) for C2 = 2. Disordered initial 
conditions: (a): C4 = 0.06; (b): C4 = 0.23. Ordered initial conditions: (c): C4 = 0.23; (d): C4 = 0.31; 
(e): c 4 = 0.38; (f): c 4 = 0.39; (g) c 4 = 0.57; (h): c 4 = 0.9. 
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FIG. 7. The coupling of the original time slices, indicated by solid lines, between two 'fat' 
slices which are represented by the boxes. The vertical dotted lines are the couplings, the nonzero 
elements of A~^. The coefficients B, C, ■ ■ ■ of the couplings are given in the lowest line. 
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